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ABSTRACT 

The hard X-ray (HXR) emission in solar flares comes almost exclusively from a very small part of the 
flaring region, the footpoints of magnetic loops. Using RHESSI observations of solar flare footpoints, 
we determine the radial positions and sizes of footpoints as a function of energy in six near-limb events 
to investigate the transport of flare accelerated electrons and the properties of the chromosphere. HXR 
visibility forward fitting allows to find the positions/heights and the sizes of HXR footpoints along 
and perpendicular to the magnetic field of the flaring loop at different energies in the HXR range. 
We show that in half of the analyzed events, a clear trend of decreasing height of the sources with 
energy is found. Assuming collisional thick-target transport, HXR sources are located between 600 
and 1200 km above the photosphere for photon energies between 120 and 25 keV respectively. In the 
other events, the position as a function of energy is constant within the uncertainties. The vertical 
sizes (along the path of electron propagation) range from 1.3 to 8 arcseconds which is up to a factor 
4 larger than predicted by the thick-target model even in events where the positions/heights of HXR 
sources are consistent with the collisional thick-target model. Magnetic mirroring, collisional pitch 
angle scattering and X-ray albedo are discussed as potential explanations of the findings. 
Subject headings: Sun: flares - Sun: X-rays, 7-rays - Sun: Chromosphere - Acceleration of particles 



1. INTRODUCTION 

In the traditional flare model particles are accelerated 
in the corona then precipitate along the field lines of 
a magnetic loop to the chromosphere where they are 
stopped producing bremsstrahlung emissi on in the pro - 
cess. In the classical thick-target model (|Brown|[l971[ ). 
it is assumed that collisional interaction of fast electrons 
with the ambient plasma leads to energy-loss, while other 
mechanisms such as pitch angle scattering or mirroring of 
the electrons in a converging magnetic field are neglected. 
It is therefore expected that electrons with higher ener- 
gies penetrate deeper into the chromosphere before they 
are fully stopped. The stopping depth depends on the 
initial electron energy and the ambient density. Ex- 
pressed in terms of the column depth N(s) = J n(s)ds, 
where n(s) is the ambient density along the electron 
path, the stopping depth is given as N stop = Eq/2K, 
where Eq is the i nitial energy of the accelerated el ectron 
and K = 27re 4 A (IBrowiilll!)72t IBrown et al.l[2002l) . The 
Coulomb logarithm A has typical values of ~ 20 in the 
(ionized) corona and ~ 7 in the (neutral) chromosphere 
|Brownlll973l:lEmsliejll978[ ). For an electron flux distribu- 
tion F(E, s) with energy E at distance s from the point 
of injection, the observed X-ray flux at Earth is: 



12001 and more recent HXR observations using albedo 
in imaging dB attaglia ~et al.l 1201 If) and spectroscopy 
ijKontar & Brown! I2006D . For increasing density along 
s, Eq. [T]has a maximum for a given photon energy e, i.e. 
the observed HXR emission will have a maximum at a 
certain chromospheric depth, depending on energy. 

Observational evidence for height dependent HXR 
source s was found early on in stereoscopic observations 
( Kanei Il983f) and in a statis t ical w ay using Yohkoh 



n(s)A(s) 
4nR 2 



F(E,s)a(e,E)dE, (1) 



where n(s) is the density, A(s) the width of the mag- 
netic flux tube at distance s, R the Sun-Earth dis- 
tance and cr( e, E) t he an gle-averaged bremsstrahlung 
cross-section (|Haud 119971 ). The isotropic approxima- 
tion of emission is supported by statistical observations 
(jKane et alJll988t iVestrand et al.lll987tlKasparova et al.1 
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(jMatsushita et al.lll992l ). Kand (|1983l) derived that HXR 
sources for energies > 150 keV should b e at heig h ts less 
than 2500 km above the photosphere. iFletcherl (1996) 
used test particle simulations to find the expected po- 
sition as a function of energy including collisional pitch 
angle scattering and magnetic mirroring and compared 
the results with observations from Yohkoh. This study 
demonstrated the influence of a converging magnetic field 
on the height of the X-ray source, showing that the 
Yohkoh observations of 53 — 93 keV sources at average 
heights ~ 6000 km presented by Matsushita et al. (1992) 
are consistent with partial trapping of electrons in a mag- 
netic loop. However, the use of Ha flare locations as the 
reference for height est imates could be the reas on for the 

large heights found by IMatsushita et al] (11992D. 

The high spatial r esolution of RHESSI (|Lin et al.ll200l 
Harford et al.1120021) now makes it pos sible to study indi- 
vidua l events with higher accuracy ( Aschw anden et aLl 
2001 iMrozekl [200l iLiu et al l [20061) . Mor e recently 



Kontar et alT(|2008D and Kontar et al.1 (|2010l) analyzed 



a limb event using the new ly developed visibility tech- 
nique (jSchmahl et al.1 12007'). This allows for measure- 
ments of HXR source pos itions with sub-arcsecond reso - 
lution. iPrato et all ()2009f ) and iPetrosian fc Chenl (|2010f ) 
went one step further and computed the electron distri- 
butions found from inversion of the X-ray visibilities. In 
all those studies a decrease of the radial position of the 
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source s with increasing energy is found. iBrown et al.l 
(2002) showed how this can be used to determine the 
chromospheric density st ructure. This was applied by 
lAschwanden et aTl (|2002[ ) to an event on 2 002 February 
20 and to a limb event on 2004 January 6 bv lKontar et al.1 
(2010) who found their observations to be consistent 
with an exponential chromospheric density profile with 
scale height pa 150 km. The HXR source heights have 
been found to decrease from pa 1200 km at 20 keV to 
pa 700 at 160 keV. A stat istical survey of over 800 flares 
(jSaint-Hilaire et aLll2010[ ) found similar heights and that 
HXR sources appear within a relatively narrow range of 
heights of 0.5 Mm. 

Due to constantly improving analysis techniques for 
RHESSI data, it is now not only possible to deter- 
mine the position of footpoint sources with high accu- 
racy, but also the characteristic siz es at different ener - 
gies and hence at different heights (jKontar et al.l [2008). 
Moreover, it is possible to assess the extent of footpoint 
sources parallel and perpendicular t o the magnetic loop, 
as demonstrated by iKontar et alj ([2010). This opens 
up a completely new approach to the diagnostics of the 
magnetic field structure of the chromosphere and the 
physics of electron transport in the chromosphere. In 
a thick-target case, electrons with higher energies pen- 
etrate deeper into the chromosphere. In a converging 
magnetic loop, the higher energetic electrons will there- 
fore penetrate into regions with higher magnetic field 
strength. This should be reflected in the horizontal 
size of the X-ray source which is expected to decrease 
with energy. Such a behavior was indee d observed in 
the event analyzed bv lKontar et al.l (|2008t ). who also in- 
ferred a magnetic scale height of the order of 300 km. 
At the same time the vertical extent of a source should 
be determined by the same p hysics of electron trans- 
port. However, IKontar et al.l (j2010) find that the ver- 
tical sizes of HXR sources are inconsistent with the sizes 
expected from the collisional thick-target. While they 
interpret this in terms of a multi-threaded chromosphere 
with different density profiles along the different threads, 
alternative explanations within a monolithic chromo- 
sphere framework such as collisional pitch angle scatter- 
ing or magnetic mirroring are possible. Another effect 
that influences the observed positions and sizes of X-ray 
source s is X-ray albedo , as sh own by IKontar fc Jeffrey! 
(|2010n . iBattaglia et al.l (|2011h went a step further and 
demonstrated how the measured full width half maxi- 
mum (FWHM) in simulated and observed maps can be 
used to constrain the true source size and the directivity. 

In this paper, we investigate the energy dependence 
of not only the position, but the sizes of footpoints in a 
carefully selected sample of six limb events observed by 
RHESSI. We used visibility forward fitting to find the 
positions and FWHM sizes of footpoints as a function of 
energy. We show that in half of the observed sources, the 
position decreases with energy and can be interpreted in 
a simple thick-target approach. In the other events the 
position is constant within the uncertainties. Contrary 
to the predictions from the collisional thick-target model, 
the vertical sizes along the electron path are a factor of 
2-4 larger and weakly decrease with energy for all flares 
analyzed. Possible explanations include collisional pitch 
angle scattering, magnetic mirroring and X-ray albedo. 



2. FLARE SELECTION 

The primary selection criterion were events with unam- 
biguous footpoint sources observed to high energies. We 
restricted our selection to events of GOES class M and 
above in order to have high enough count rates. At the 
same time events with strong pulse pile-up (live -time less 
than « 85%) were excluded (| Smith et al JI2002! ) . Pile-up 
appears when two or more low energy photons, mostly 
from the coronal source, are detected as a single pho- 
ton with higher energy. It is of less concern in footpoint 
analysis, as the pile-up signature in images appears as ar- 
tificial high energy source in the corona where the bulk of 
the soft X-ray emission originates. Except for the 2005 
August 22 the attenuator state during the events pre- 
sented here was 1. Thus the peak of the pile- up is at 24 
keV, below our analyzed energies. In the 2005 August 22 
(attenuator state 3), the live-time was high enough for 
pile-up to be minimal. To test the effect potential pile-up 
counts would have on the visibility forward fit, we fitted 
the two footpoints in the 2005 July 13 event in an energy- 
band of 16-40 keV in which emission from the coronal 
source is present, mimicking pile-up. It is found that the 
positions are shifted by less than 0.5 arcsecond and the 
FWHM change by less than 0.2 arcsec, both effects being 
smaller than the uncertainties in the fit parameters even 
though the intensity of the coronal emission per unit area 
was only a factor of 2 less than the footpoint emission 
per area. Further, to reduce projection effects and find 
events with a clear morphology, the event selection was 
restricted to limb events (angular offset from Sun center 
larger than 700 arcsec). This also minimizes the influ- 
ence of X-ray albedo on the source positions and sizes 
(Section 03). Quick look images prov ided by the HESSI 
Expe rimental Data Center (HEDC0, iSaint-Hilaire et all 
2002) were used to search for events with one or two foot- 
points visible up to at least about 80 keV. Events with 
more than two footpoints or o ther complex stru ctures 
such as the 2002 July 23 event (|Emslie et al.l 12003) were 
excluded. Finally, six events satisfying aforementioned 
criteria were selected (The key parameters of the flares 
are listed in Table CQ). 

Figure 00 gives the temporal and morphological 
overview of the events. Contours at 50 %, 70% and 90 % 
of the maximum emission in CLEAN images are shown 
for energy bands 6-12 keV and 35-55 keV. The 6-12 keV 
emission represents the thermal (coronal) emission, while 
the 35-55 keV emission indicates the footpoint emission. 
Two of the events only have one footpoint. Two events 
have one strong and one weak footpoint while in the other 
two events, the intensities of HXR footpoint emission are 
comparable. The lightcurves of corrected count rates in 
the 6-12 keV and 25-50 keV energy bands are presented 
underneath the images. 

3. DATA ANALYSIS 

For the analysis we chose a time interval of 1 minute 
during which the footpoint emission is observed to high- 
est energies in the quicklook images on HEDC (except in 
the 2002 April 26 event during which attenuator state 
changes only allowed an interval as long as 12 sec). 
This does not necessarily coincide with the peak of the 

1 Unavailable now as it went off-line at the beginning of 2010 
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Figure 1. Images and lightcurves of the selected events. Contours of 50%, 70% and 90% of the maximum emission in CLEAN images 
using detectors 3-8 (1-8 for the high energies in the 26 April event) are given for energy ranges 6-12 keV (red, thermal emission) and 
35-55 keV (blue, non-thermal emission). Lightcurves of corrected count rates in the 6-12 keV (red) and 25-50 keV (blue) energy bands are 
displayed underneath each image. The dashed lines mark the analyzed time interval. 



Table 1 

Dates, times, GOES class, electron spectral index <5 and cosine of 
heliocentric angle 9 of the flaring region for the analyzed flares. The times 
given are the start times of the 1-minute observation interval (Section (3j, 
except for the 2003 April 26 event for which the full interval (less than 1 
minute) is given. 



Date Start time [UT] GOES class spectral index 8 cos 9 



2002 May 31 00:06:42 M2.5 3.5 0.04 

2003 Apr 26 23:39:28-23:39:40 Ml.l 3.6 0.16 

2004 Jan 7 10:21:20 M8.3 3.4 0.28 

2005 Jul 13 14:15:00 M5.1 4.1 0.2 
2005 Aug 22 17:09:46 M6.3 3.7 0.47 
2005 Aug 23 14:37:00 M3.0 4.5 0.29 
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HXR emission. Shorter time intervals generally result in 
poorer signal-to-noise ratio for HXR visibilities. Further, 
the longer time interval improves the aspect phase cov- 
erage because of the precession of the RHESSI spin axis 
(sec 13.11) a nd also helps to minimize the effect of data 
gaps (Hur ford et al.ll2002t ). On the other hand it is often 
observed that footp oints move spatially over the course 
of on e minute (e.g. jGrigis & Bend [2005; Fle tcher et al.l 
2004; IKrucker et aU l2003h which might affect the ob- 
served positions and sizes. In the events presented here, 
the motion of the sources is found to be less than one 
arcsec during the observed time interval. This is of the 
order of uncertainties in the measurements of the posi- 
tions and sizes. Therefore, this effect can be neglected. 

Images at several energy bands in the non-thermal en- 
ergy range were made to find the positions and the sizes 
of the footpoints. The energy-binning was chosen to pro- 
vide 5-6 images in the hard X-ray domain per event, 
star ting from 25 keV (Detector 2 sensitivity threshold, 
see iSmith et al.l l2002f) up to the highest observed en- 
ergies (100 - 200 keV, typically). A finer energy bin- 
ning would be desirable but would lead to insufficient 
number of cou nts in the individual energy bands . We 
used C LEAN (iHogboml 119741: iHurford et al.l 12002ft and 
Pixon (|Pina & Puetterl Il992t iMetcalf et aT[ [1996) for a 
first impression of the sources, their intensity and the 
shape. The positions and sizes as a function of energy 
were then f ound using the technique of visibility forward 
fitting (e g.lBattaglia et alJl201lHDennis fc Pernakl l2009t 
IXu et al.ll2~008t Hannah et al.ll2008ft~ 

3.1. Visibility forward fitting 

Visibilities are a concept widely used in radio astron- 
omy. In rece nt years it has been adapted to be use d 
with RHESSI (IHurford et al.ll200l ISchmahl et al.ll2007ft . 
HXR visibility forward fitting is ideally suited to find po- 
sitions and sizes of flare sources for several reasons. X-ray 
visibilities V(u, v) are the 2D spatial Fourier components 
of the X-ray distribution I(x, y): 

V(u, v) = / / I(x,y) exp[2ni(ux + vy)]dxdy, (2) 

J x J y 

so the total flux, position and size of an X-ray source 
(comp. Eq. [1]) can be straightforwardly expressed 
through the 0th, 1st and 2nd moment of the X-ray dis- 
tribution. Measurements of the 1st and 2nd moment 
from reconstructed images (e.g. CLEAN, Pixon) are 
impractical since the result is dependent on the selec- 
tion of an image region. Because RHESSI images suffer 
from reconstruction noise and are hampered by imag- 
ing artefacts such as side lobes from the CLEAN beam, 
so will the moments. Visibility forward fitting does not 
require the reconstruction of an image itself and is there- 
fore th e most direct way to find a measure of the mo- 
ments (iBattaglia et al.ll2011[ ). If the observed source has 
a Gaussian shape then the forward fit parameters will 
represent the true moments. Further, statistical errors 
for the fit parameters are computed from the visibility 
errors making visibility forward fitting the only method 
that provides error estimates for the measured parame- 
ters. Figure [5] gives an example of a visibility forward 
fit. All detectors were used since we only analyze ener- 
gies above the sensitivity threshold of detector 2 (> 25 



keV). On the left-hand side a Pixon map is shown, over- 
laid with the 50%, 70% and 90 % contours of the source 
model (two circular Gaussians). The top panel on the 
right-hand side shows the observed visibility amplitudes 
with statistical uncertainties and the fitted model. The 
bottom panel on the right-hand side displays the normal- 
ized residuals of the visibility amplitudes. Note that the 
model is fitted using all V(u, v), not only the amplitudes 
\V(u, v) \ shown in Figure[2] Apart from the source shape 
given by the model, the number of roll-bins i.e. the spa- 
tial coverage for which visibilities are computed, has to 
be chosen. A large number of roll bins will provide finer 
spatial coverage. If used, especially for the finer grids 
(detectors 1-3), potential ellipticity of the source can be 
better assessed and fitted, so the visibilities might even 
reflect smaller source structures. At the same time, the 
errors of the individual visibilities tend to be larger for 
small count rates as the same number of counts is spread 
over a larger number of visibilities. A smaller number 
of roll bins provides coarser spatial coverage but with 
smaller errors of the individual visibilities. However, the 
choice of roll bin number does not affect the inferred 
forward fit parameters (flux, position and size), as tests 
using different settings proved. 

3.2. Spectroscopy 

Full Sun spectra were analyzed for all events to gather 
information about the thermal and non-thermal contri- 
bution to the emission. This allows to determine the 
lower limit of non-thermal energies that can be analyzed 
in imaging without the risk of contamination from the 
coronal source. The spectra of all events were fitted with 
a thermal model plus a thick-target power-law model 
(Fig. [3|). This provides the electron spectral index S 
of accelerated or into the thick-target injected electrons 
that is used as a parameter in fitting the chromospheric 
density. The comparison of the flux from full-Sun spec- 
tra with the total flux from the footpoints as found in 
visibility forward fitting is an independent test of the 
consistency of the analysis (Fig. [3]). Assuming the HXR 
emission comes predominately from the footpoints, the 
total flux from the individual sources equals the flux from 
the full-Sun spectra if it has been accounted for correctly. 
The full-Sun spectra are shown in Fig. [3] with the flux 
from visibility forward fitting overlaid. Generally there 
is a good agreement between the full Sun flux and the to- 
tal footpoint flux found in visibility forward fitting. The 
exceptions are the points at the lowest energies in the 
events of 2005 July 13 and 2005 August 23, which could 
be related to some non-footpoint emission from the coro- 
nal source, albedo or the loop itself. 

4. RESULTS 

Although all flares have a number of common features, 
some aspects are highly individual and it is difficult to 
make general statements that hold for all. We there- 
fore make some general remarks on the results and ex- 
plain how the chromospheric density was inferred from 
the measurement of the position as a function of energy. 
Then, each event is discussed individually in more de- 
tail. Figure U presents the measurements of radial posi- 
tion and source size for all events. The first and third 
column illustrate the measured radial distance (2D dis- 
tance measured from the solar disk centre) as a function 



Hard X-ray footpoint sizes and positions 



5 



(D 

















i 




ft 




k 






Figure 2. Illustration of visibility forward fit. Left: Pixon map overlaid with the 50%, 70% and 90% contours of the source model (two 
circular Gaussians). Top right: Measured visibility amplitudes (black dots) with statistical uncertainties (blue) as a function of position 
angle for each RHESSI detector (1 to 9). Bottom right: Normalized residuals of the visibility amplitudes. 
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Figure 3. Full Sun spectra, fitted with a thermal model (green) and a non-thermal power-law (purple). The blue lines indicate the flux 
from the individual sources found in visibility forward fitting. The red lines give the total of the individual sources in the case of two 
footpoints. 



of energy. The second and fourth column show the size 
as a function of energy. 

4.1. Radial position and density 

Visibility forward fitting returns the x — y position (1st 
moment) of the fitted source including errors as a func- 
tion of energy (x(e),y(e)). The radial position can read- 
ily be found as r(e) = \J x{e) 2 + y(e) 2 . In the classical 
thick-target model HXR emission at higher energies is 
expected to originate deeper in the chromosphere. This 
can be used to find the chromospheric density profile by 
fitting a model to the observed positions. We use the 
same approach and general assumptions as lKontar et al.l 



(201(j). A hydrostatic density profile of the form 

h 



n(h) = n c + no cxp [ — — 
ho 



(3) 



is assumed where ho is the density scale height, no the 
photospheric density and n c = 9 x 10 9 cm -3 is the coro- 
nal density, taken as constant. The density profile is 
illustrated in Fig. [5] Based on iVernazza et al.l (|1981[ ) 
we chose the fixed value of no = 1.16 x 10 17 cm -3 for 
h = (photospheric level). Assuming a vertical loop, 
the height of the source above the photosphere, h(e) can 
be expressed through the radial distance from the Sun 
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Figure 4. Results from visibility forward fitting for all events. First and third column: Radial distance as a function of energy. Second 
and fourth column: FWHM as a function of energy. 



center as h(e) = r(e) — r$ where r is the radial distance 
that corresponds to the photospheric height h = 0. The 
parameters kg and r$ can be found by forward fitting the 
measured radial distances r(e) with the density model in 
Eq.El 

4.2. Size and shape 

The second parameter of interest is the source size. In 
visibility forward fitting one or two Gaussian sources 



I(x, y) ~ exp 



x 

2^ 



y 

2al 



(4) 



are fitted, where a a and <7& are the standard devia- 
tions. Visibility forward fit returns the FWHM and 
the eccentricity e from which the major and minor axes 
(a = 2V21n2cr a and b = 2V2hi 2<7f,) as a function of 
energy can be found: a(e) = FWHM(e)(l — e(e) 2 )~ 1//4 
and 6(e) = FWHM(e)(l - e(e) 2 ) 1/4 . In the case of a 



circular Gaussian fit, the FWHM is simply the diam- 
eter at half the peak-flux level and a = b. If the fitted 
source is a true Gaussian in reality, the size is directly re- 
lated to the second mo ment of the X-ray flux distribution 
(jBattagha et al.ll20llD . Depending on the loop geometry 
the major and minor axes are the sizes perpendicular 
and parallel to th e magnetic field of t he loop (as in the 
event analyzed bv lKontar et al.ll2010f ). In this case, the 
minor axis is measured radially from Sun center and rep- 
resents the size vertical to the photosphere while the ma- 
jor axis is measured perpendicular to the radial direction, 
i.e. parallel to the photosphere. In a simple thick-target 
model, the vertical size depends on the density structure 
and the electron spectral index. Using the density model 
found as described in Section I4.ll the expected X-ray 
source profile along the loop can be calculated for differ- 
ent energies. This is illustrated in Fig. [5] The expected 
photon flux as a function of height above the photosphere 
is shown for different energy ranges. The vertical extent 
of the source is given as the width of the curve at the level 
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of half the peak-flux. It can also be seen how the position 
of the maximum is shifted to lower heights for increasing 
energy. Fig. [6] illustrates the expected vertical extend in 




Figure 5. Illustration of the photon flux as a function of height 
above the chromosphere for different energy bands (explained in the 
legend) in a simple thick-target model. The width of the curves at 
the level of half the peak flux gives the vertical size of the source. 
Inset: Density profile. 

the thick-target model compared to the observations for 
selected events. In the following we are discussing the 
individual events in more detail. 

4.3. 2002 May 31 

This is a limb event with only one footpoint (Fig. [TJ 
top left). The event was best fitted with a single ellip- 
tical Gaussian. A decrease of the radial position with 
energy is observed, although the individual errors on the 
position are rather large (Fig. 0|). We fitted the density 
model described in Eq. [3] finding a density scale height 
of ho = 204 ± 37 km. The observations also suggest a 
decreasing FWHM with energy, although with large un- 
certainties. This is the only event that had a clearly ellip- 
tical footpoint shape. The major and minor axes as well 



T 



30 40 50 60 70 H 
Energy [kev] 

26/04/2003 north 



20 30 



Energy [keV] 
26/04/2003 south 



40 50 60 70 H 
Energy [kev] 

13/07/2005 s 



+- 









Figure 6. Vertical extent (FWHM in case of circular Gaussian 
sources, FWHM of minor axis in the case of the 2002 May 31) 
as a function of energy. Th e red points indicate the projection 
corrected extents (Sect. 15.111 . The blue line gives the expected 
vertical extent from the density model that was found by fitting 
the radial positions of the respective events. 



as the eccentricity are displayed in Fig. [7] The major 
axis is oriented along the limb. The morphology suggests 
that the major axis of the ellipse is perpendicular to the 
magnetic field of the loop, while the minor axis is parallel 
to the field. The extent of the major axis decreases with 
increasing energy, which can be interpreted as the signa- 
ture of electron transport in a converging magnetic field. 
The minor axis is constant as a function of energy and 
much larger than the size expected from the thick-target 
model. 
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Figure 7. Major and minor axes (top), eccentricity (bottom) in 
the event of 2002 May 31 which was fitted with an ellipse. 



4.4. 2003 April 26 

This event has two very compact, close footpoint 
sources (named north and south in Fig. 2]) that were 
best fitted with two circular Gaussian sources. In both 
sources, a tendency to a smaller radial distance with in- 
creasing energy can be observed. The fitted density scale 
heights arc ho = 132 ± 38 km for the northern foot- 
point and ho — 181 ± 40 km for the southern footpoint. 
The FWHM of the northern footpoint is constant as a 
function of energy. The size of the southern footpoint 
decreases as a function of energy except for the mea- 
surement at highest energies, which might be due to low 
count rate in the source. Interestingly, the size at 40 keV 
to 60 keV of the northern footpoint is consistent with 
the thick-target prediction, but not at the other energies. 
This indicates that geometrical effects such as projection 
or footpoint motion are negligible in this case since they 
would affect all energy ranges equally. 

4.5. 2004 January 7 

This flare happene d one day aft e r the 2004 January 
6 event analyzed by iKontar et all 1)20081 l2010f ) in the 
same active region. It has two footpoints although one 
is much fainter than the other. A model of two circular 
Gaussians is fitted to include the emission of the weak 
footpoint but only the results for the stronger footpoint 
are shown in Fig. U (top right) due to the large uncertain- 
ties resulting from the small count rates in the weaker 
footpoint. The density scale height was found to be 
ha = 230 ± 40 km. The size is constant as a function 
of energy except in the lowest energy band in which the 
size is considerably larger. Although the spectrum sug- 
gests purely non-thermal emission in this energy band, 
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it cannot be ruled out entirely that there is still some 
coronal emission which might affect the fitted size. 

4.6. 2005 July 13 

This event has two footpoints (named north and south 
in Fig. [4} of equal intensity that are best fitted with two 
circular Gaussians. The radial position of the southern 
footpoint clearly decreases with energy. A density scale 
height of ho = 235 ± 56 km is found in this case. The 
height of the northern footpoint source also decreases as 
a function of energy, except in the lowest energy band. 
The size of both footpoints can be considered constant. 

4.7. 2005 August 22 

The event has two footpoints with the northern foot- 
point being much fainter than the southern footpoint. 
Two circular Gaussians were fitted, but only the results 
for the southern footpoint are shown in Fig. 2] There 
is an indication for a decrease of position with energy 
followed by an increase, although the uncertainties are 
rather large. FWHM size is approximately constant. 

4.8. 2005 August 23 

Only one distinct footpoint is measurable in this event 
although imaging over a larger energy band suggests the 
existence of a second, very faint footpoint. One circular 
Gaussian source was the best model for this case. Within 
the uncertainties, radial position as well as FWHM are 
constant. 

5. PHYSICAL INTERPRETATION OF THE OBSERVATIONS 

While some of the described events display a clear de- 
crease of radial distance with energy (similar to previous 
results) constant positions are also observed. In the clas- 
sical thick-target model, one would expect a decrease of 
the height of a source with increasing energy. There are 
a number of effects that can influence the measured po- 
sitions and sizes. 

5.1. Projection effects 

The most obvious is projection effects. In observa- 
tions we measure the position of the source in radial di- 
rection from the center of the Sun. For sources at or 
close to the limb, these positions are directly related to 
the height of the source above the photosphere. Closer 
to the center of the Sun, sources at different heights in 
the chromosphere will be seen in projection on top of 
each-other. This would result in an observed constant 
radial position as a function of energy. To avoid this 
we focused on near limb events. However, even for near 
limb events, an effect due to the heliocentric angle is 
expected. Table 1 lists the cosine of the heliocentric an- 
gles of the flaring region. The 2005 August 22 event is 
the furthest away from the limb at a heliocentric angle 
of = 62°, corresponding to [i = cos(#) = 0.47. This 

introduces a factor 1 / y/l — fj, 2 to the effective heights 

htme = h fated/ \/i ~~ A* 2 which could affect the density 
fits. While no density fit was possible for the 2005 Au- 
gust 22 event (y/l — fjp — 0.88) because the positions as 
a function of energy were constant within uncertainty, 
in all the other events \/l — /i 2 > 0.96 and the result- 
ing effect is smaller than the uncertainties of the fits. 



Therefore, the influence of the heliocentric angle on the 
measured source heights and derived densities is negligi- 
ble. The contribution of projection to the measured sizes 
can be estimated using a simple geometrical model. For 
a source with fitted horizontal extent (major axis) a the 
fitted vertical extent (minor axis) b is composed of: 

b = b^ue x M 2 ) +ax fi (5) 

where btme is the true vertical extent. For events exactly 
at the limb /i = 0, thus b = btme- Assuming a circular 
shape of the footpoint at a given height the extent per- 
pendicular to the radial direction, which is unaffected 
by projection effects, can be used to estimate the true 
vertical size. Figure |6] illustrates the effect. 

In all of the observed events, the measured source size 
is larger by at least a factor of three compared to the 
expected size in a thick-target model (Fig. [6]). Even 
when projection effects are included, the sizes are still 
more than a factor of 2 larger than expected from the 
thick target model. Currently there is only one expla- 
nation i.e. multi-threaded loop density structure that 
has been investigated in more detail and that was used 
to explain the obse rved size in the event analyzed by 
iKontar et al.l (|2010l ) . Another possibility related to the 
density structure is a double-exponential density struc- 
ture with a sec ond, larger scale-height at higher altitudes 
as proposed bv lSaint-Hilaire et al.l (|2010[) . This might af- 
fect the height of the sources as a function of energy and 
possibly the sizes. However, it is unlikely that the ef- 
fect will be large enough to explain the observed sizes. 
Moreover, in terms of the positions in the events pre- 
sented here a density structure with a single scale-height 
is sufficient to explain the observed function of radial dis- 
tance versus energy. However there are other effects such 
as magnetic mirroring, pitch angle scattering or X-ray 
albedo that are expected to affect the size of the source. 

5.2. Magnetic mirroring 

In the standard thick-target model, magnetic mirror- 
ing is neglected. However, in a converging magnetic field 
it is expected that some particles are mirrored back from 
the footpoints, if they were injected at a pitch angle rel- 
ative to the magnetic field lines. The mirroring point 
depends on the initial pitch angle of the electrons but 
is independent on the electron energy. The bulk of the 
emission is expected to originate from the densest part 
of the chromosphere to which the electrons are able to 
penetrate. A single mirroring point would lead to a con- 
stant position as a function of energy such as observed in 
the 2005 August 23 event. A behavior as in the stronger 
source of the 2003 April 26 event could also be envis- 
aged. If the stopping depth for low energetic electrons 
is substantially higher than the mirroring point, those 
electrons will encounter a thick-target, while higher en- 
ergetic electrons will be mirrored. Therefore, a decrease 
of radial position with energy will be observed at low en- 
ergies, a constant position at higher energies. The size of 
HXR sources will depend on the pitch angle spread of the 
electrons. Field aligned electrons will penetrate deeper 
while electrons with a large pitch angle will be mirrored 
higher. This would make the source larger. On the other 
hand, in the extreme case of injection at a 90° angle to 
the magnetic field the electrons will stay at the height of 
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the injection point, gradually losing energy. This would 
result in a constant position and a source size determined 
by the size of the acceleration region. 

5.3. Collisional pitch angle scattering 

Another effect that might not be negligible is collisional 
pitch angle sc attering. It is likely to increase the size 
of the source (|Conwayl 120001 ) . but this increase is not 
expected to be large enough to explain the observations. 
Additional scatteri ng due to various plasma waves can 
be anticipated (e.g.lBian et al.ll2010l : lHannah et al.ll200li 
iStepanov fe Tsap|l2002D ." While this can be substantial, it 
is difficult to quantify the level of turbulence in a flaring 
atmosphere. 

5.4. X-ray albedo 

Photon backscattering (albedo) from the photosphere 
can change the observe d po sition and size as sho wn by 
iKontar fe Jeffrevi (pOlOh and lBattaglia et all (pOlTI ) . The 
additional albedo flux results in a shift of the position 
radially toward the disk center and in an increase of 
the observed source size by several arcseconds. Since 
albedo is energy dependent (with the strongest contribu- 
tion between around 30 and 50 keV), the effect on the 
position and the size will be most pronounced at those 
inter mediate energies and sm aller at low and high ener- 
gies. IKontar Sz JeffrevI ([2010D further find that this effect 
can be important even for large heliocentric angles, i.e. 
events near the limb. Being the closest to the solar disk 
(at a cosine of the heliocentric angle of fj, = 0.47), the 
observations of the event of 2005 August 22 might be 
affected by albedo. The positions would be shifted radi- 
ally toward the disk center, most pronounced at energies 
between 30 and 50 keV. This could explain the observed 
pattern (Fig. 21 bottom left). 

6. DISCUSSION AND CONCLUSIONS 

In this work we performed the first multi-event study 
of both source positions and sizes of flare footpoints as 
a function of energy. The radial positions as a function 
of energy in 4 out of 6 events follows a decreasing trend 
with increasing energy as expected in the classical thick- 
target model. The positions were fitted with an expo- 
nential density model, finding scale heights between 132 
and 235 km, consistent with previous studies. The other 
two events show a weak dependence of the positions on 
energy. The vertical sizes range from 1.3 up to 8 arcsec, 
therefore being up to a factor 3 larger than expected from 
the simple thick-target model, even for the events with 
clear dependence of the position on the energy and after 
correcting for projection effects due to the heliocentric 
angle (see Fig. |6|). While it has been argued, based on 
traditional imaging techniques, that R HESSI might sim- 
ply u nder-resolve footpoint sizes (e.g. iDennis fc Pernakl 
2009) , other studies suggest that X-ray f ootpoint sizes of 
the order of several arcse conds are real (iBattagli a et al.l 
120111: IKontar et al.ll20To| ). This is supported by negli- 
gible modulation in the finest RHESSI grid (grid one), 
which gives a lower limit of the source size of 3.9 arcsec 
in the events presented here. Currently, the only expla- 
nation for this in the context of a simple thick-target 
model is a multi-threaded loop density structure. Here 
we discuss other mechanisms which might explain this 



intriguing finding. In the thick-target model, electrons 
are injected parallel to the magnetic field lines. In real- 
ity, injection at a pitch angle relative to the field lines can 
be expected. In a converging magnetic field, the initial 
pitch angle distribution will be modified and magnetic 
mirroring will play a non-negligible role. Depending on 
the height of the mirroring point, the observed height 
as a function of energy will be constant. This can ex- 
plain the observations of the position in two events. At 
the same time, the size is expected to be affected, as 
well. In addition collisional pitch angle scattering, as 
well as wave-particle interactions will further modify the 
pitch angle and therefore the observed size, although col- 
lisional pitch angle scattering alone will not be sufficient 
to enhance the size. A third possibility is the influence 
of X-ray albedo. Although the effect is expected to be 
small in limb events, it might affect the measurements 
in the event of 22 August, which is closest to the disk 
center. 

While the presented events are very individual, a gen- 
eral trend can still be found when averaging the sizes. 
Figure [8] shows the average vertical and horizontal size 
as a function of energy. In the case of a circular Gaus- 
sian source, the vertical and horizontal sizes are identical 
and equal to the FWHM, while in the case of the 2002 
May 31 event, the minor and major axis were used. The 
figure suggests a decrease of both dimensions with en- 
ergy. In order to understand the physics and explain 
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Figure 8. Average horizontal (top) and vertical (bottom, projec- 
tion corrected) extent of the sources. The error bars indicate the 
standard deviation of the scatter of the individual points. 

the observations, one needs to carefully address both the 
position and the sizes since most of the effects discussed 
are expected to affect both, size and position. Test parti- 
cle simulations including magnetic field, collisional pitch 
angle scattering as well as other effects such as pitch an- 
gle scattering due to turbulence are necessary and have 
to be compared with the observations. Such a complete 
treatment will make it possible to use X-ray footpoints as 
an independent means of determining the chromospheric 
magnetic field and density structure. 
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